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We present a near-infrared extinction map of a large region (approximately 2 200deg ) covering the Orion, the 
Monoceros R2, the Rosette, and the Canis Major molecular clouds. We used robust and optimal methods to map 
the dust column density in the near-infrared (Nicer and Nicest) towards ~ 19 million stars of the Two Micron All 
Sky Survey (2MASS) point source catalog. Over the relevant regions of the field, we reached a 1-a error of 0.03 mag in the 
K-hand extinction with a resolution of 3 arcmin. We measured the cloud distances by comparing the observed density of 
foreground stars with the prediction of galactic models, thus obtaining dorionA = (371 ±10) pc, dorionB ~ (398±12) pc, 
dMonR2 = (905 ± 37) pc, dRosottc = (1330 ± 48) pc, and dcMa ~ (1150 ± 64) pc, values that compare very well with 
independent estimates. 

Key words. ISM: clouds, dust, extinction, ISM: structure, ISM: individual objects: Orion molecular complex, ISM: 
individual objects: Mon R2 molecular complex. Methods: data analysis 



1. Introduction 

In a series of papers, we have applied an optimized 
multi-band techniq ue dubbed Near-Infra r ed Co lor Excess 
Revisited (Nicer Lombard! & Alves 2001 hereafter 
Paper 0) to measure dust extinction and investigate the 
structure of nearby molecular da rk clouds using the Tw o 



Micron Ah Sky Survey (2MASS; Kleinmann et al.||1994 ). 
The main aim of our coordinated study is to investigate 
the large-scale structure of these objects and to clarify the 
link between the global physical properties of molecular 
clouds and their ability to form stars. Previously , we con- 
sidered the Pipe nebula (see Lombardi et al. 2006 hereafter 
Paper I), th e Ophiuchus and Lupus complexes (Lombardi 



et al. 



2008 



hereafter Pap er II), and the Taurus , Perseus 

and Cialifornia complexes ( Lombardi et al.| 2010 hereafter 
Paper III). We now present an analysis oi a large region 
covering more than 2 200 square degrees, centered around 
Orion. This region includes Orion A and B, A Orionis, Mon 
R2, Rosette, and Canis Major. An overview of a subset of 
the wide field extinction map presented in this paper, su- 
perimposed on an optical image of the sky, is presented in 
Fig.[T} 

Near-infrared dust extinction measurement techniques 
present several advantages w ith respect to other co lumn 
density tracers. As shown by Goodman et al. (20091), ob- 
servations of dust are a better column density tracer than 
observations of molecular gas (CO), and observations of 
dust extinction in particular provide more robust measure- 
ments of column density than observations of dust emission, 
mainly because of the dependence of the latter measure- 
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ments on uncertain knowledge of dust temperatures and 
emissivities. Additionally, the sensitivity reached by near- 
infrared dust extinction techniques, and by Nicer in par- 
ticular, is essential to investigate the low-density regions of 
molecular clouds (which are often below the column density 
threshold required for the detection of the CO molecule) 
and therefore to estimate the mass of the diffuse gas that 
acts as a pressure boundary around the clumps. Similar 
to Paper III, we use for some key ana lyses in this paper 
the improved Nicest method (Lombardi 2009), designed 
to cope better with the unresolved inhomogeneities present 
in the high-column density regions of the maps. 

Orion is probably the best studied molecular cloud in 
the sky. The complex comprises Hii and Hi regions super- 
imposed on colder, massive H2 clouds with active formation 
of both low and high mass stars. The area studied here 
contains the Orion OB association, which is split in sev- 
eral subgroups with ages from 2 to 12 Myr. It is estimated 
that in the last 12 Myr there have been 10 to 20 supernova 
explosions (Bally 2008) that shaped the gas in the region. 



More than a century ago, Barnard discovered a large arc of 
Hjj emission around the eastern part of Orion (see Fig. fl]) . 
More recently, Barnard's loop has been linked to Eridanus 
loop, and the two have been identified as a superbubble that 
is expandin g into denser regions with a mean velocity of 10- 
20 km s"^ ( Mac Low fc McCray|l988 ). Orion includes two 
giant molecular molecular clouds, Orion A and Orion B, 
the spectacular A Orionis bubble, and a large number of 
smaller cometary clouds. Both Orion A and B are observed 
in projection inside Barnard's loop and most likely have 
been shaped by the supernova explosions, stellar winds, and 
Hii regions generated by OB stars in the area; for exam- 
ple, note the shell originating from the south-east extrem- 
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Fig. 1. An optical image of Orion (Wei-Hao Wang, IfA, University of Hawaii) with the extinction map presented in this 
paper superimposed in green. The complementarity between the red Hii regions and the green H2 ones is evident. Note 



that the online figure is interactive: extinction and labels can be toggled by clicking on the respective boxes when 
using Adobe® Acrobat®. 
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ity of Orion A, around the early B star k Orionis (Saiph). 
Orion A and B host rich clusters of young (^ 2 Myr) stars, 
the Orion Nebula Cluster (ONC), NGC 2024 (also known 
as the Flame Nebula), NGC 2071, and NGC 2068. Just 
southwest of NGC 2024 is the famous Horsehead Nebula 
clearly visible as a protrusion of the green extinction map 
in Fig. [1] 

Around the "head" of Orion, the 08 star A Orionis, 
there is a ring of dark clouds, 60 pc diameter, known for 
almost a century. The region hosts several young stars (in- 
cluding 11 OB stars very close to A Orionis), and their 
spatial and age distributions show that originally star for- 
mation occurred in an elongated giant molecular cloud 



( Mathieu|2008[|Duerr et al.|1982[ ). Most likely, the explosion 
of a supernova coupled with stellar winds and Hii regions 
destroyed the dense central core, created an ionized bubble 
and a molecular shell visible as a ring. Star formation still 
continues in remnant dark clouds distant from the original 
core. 

The Monoceros R2 region is distinguished by a chain 
of reflection nebulae that extend over 2° on the sky. The 
nomenclature "Mon R2" indicates the second association 
of reflection n ebula in the constellation Monoceros Ivan denl 
Bergh (1966). The region is also well known for its dark 
nebula, whicn is clearly visible on top of the reflection neb- 
ulae and the field stars. Clusters of newly formed stars are 
present in its core and in the GGD 12-15 region; smaller 
cores are also present in the field. The cloud is estimated to 
be at a dista nce significantly higher than the Orion nebula, 
830 ± 50 pc (iHerbst & Racine||1976 1 . 



Rosette is known to amateur astronomers as one of the 
most spectacular nebulae in the sky. The region is charac- 
terized by an expanding Hll region interacting with a giant 
molecular cloud. The Hii region has been photodissociated 
by NGC 2244, a cluster containing more than 30 high-mass 
OB stars. Judging from the age of NGC 2244, the Rosette 
molecular cloud has been actively forming stars for the last 
2 to 3 Myr, and most likely will continue doing so, as there is 
still sufficient molec ular material placed in a he avily stimu- 
lated environment ( Roman- Ziifiiga & Ladar2008). The dis- 
tance to NGC 2244, and therefore to the molecular cloud, is 
still controversial: recen t estimates ra nge between 1 390 p c 
(Hensbergeet al.||2000|) and 1 670 pc (iPark & Sung"2002). 



The star-forming region in Canis Major is characterized 
by a concentrated group of early type stars (forming the 
CMa OBl/Rl associations) and by an arc-shaped molec- 
ular cloud (probably produced by a supernova explosion). 
Similarly to Rosette, Canis Major presents an interface be- 
tween the Hii region and the neutral gas which shows up 
as a north-south oriented ridge. Recent distance determi- 
nations of C anis Major seem to con verge around a distance 
of 1 000 pc ( [Gregorio He"te^[2008| . 

This paper is organized as follows. In Sect. [2] we briefly 
describe the technique used to map the dust and we present 
the main results obtained. Section [3] is devoted to an in- 
depth statistical analysis, and includes a discussion of the 
measured reddening laws for the various clouds, the mea- 
surements of the distances using foreground stars, the log- 
normality of the column density distributions, and the ef- 
fects of small-scale inhomogeneities. The mass estimates for 
the clouds are presented in Sect. [4] Finally, we summarize 
the results obtained in Sect.[5l 

A few plots in this paper are associated with switches 
designed to show or hide labels; these are shown as frames 
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Fig. 2. Color-color diagram of the stars in the whole field. 
The contours are logarithmically spaced, with each contour 
representing a density ten times larger than the enclos- 
ing contour; the outer contour detects single stars. Only 
stars with accurate photometry in all bands (maximum 1- 
a errors allowed 0.1 mag) have been included in this plot. 
The diagonal spread of stars is due to reddening along the 
marked reddening vector of dwarf and giant stars. The ex- 
pected colors of these stars, as predicted by.Bes sell fc Brett" 
( 1988 ) and converted int o the 2MASS ph otometric system 
using the relations from ( ;Carpenter|2001 1 , are shown in the 
41 the observed ci 



plot, and match very wel 



colors. 



in the respective captions. In order to use this feature this 
electronic document should be displayed using Adobe® 
Acrobat®. 



2. Nicer and Nicest extinction maps 

The data analysis was carried out following the technique 
presented in Paper and used also in the previous papers of 
this series, to which we refer for the details (see in particular 
Paper III). We selected reliable point sourc e detections from 
the T wo Micron All Sky Surve3p](2MASS; |Kleinmann et al 
1994) in the region 



180° < I < 240° 



-40° < & < 0° 



(1) 



This area is '-^ 2 200 square degrees and contains approxi- 
mately 19 million point sources from the 2MASS catalog. 
The region encloses many known dark molecular cloud com- 
plexes, including the Orion and Mon R2 star forming re- 
gions, the A Orionis bubble, the Rosette nebula, and the 
Canis Major complex (see Fig. l3|. 

As a preliminary step, we constructed the color-color 
diagram of the stars to check for the possible presence of 
anomalies in colors of stars. The result, shown in Fig. [2| 
displays a weak sign of bifurcation in the distribution of 



See http: //www. ipac . caltech.edu/2mass/. 
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Fig. 3. The Nicer extinction map of the Orion, Mon R2, Rosette, and Canis Major complexes. The resolution is 
FWHM = 3 arcmin. The various dashed boxes mark the regions shown in greater detail in Figs. 4l|7 Toggle labels . 
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Fig. 5. A zoom of Fig. [s] showing the A Orionis bubble, be due to the different colors of extinguished field stars and 
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Asymptotic Giant Branch (AGB) stars. However, given the 

weakness of the AGB stars contamination, we decided to 

proceed similarly to Paper III and not exclude any object. 

After the selection of a control field for the calibration of 

source colors along the reddening vector. This is likely to the intrinsic colors of stars (and their covariance matrix) we 



M. Lombardi et at: 2MASS extinction maps: IV. Orion, Mon R2, Rosette, & Canis Major 



3 -2- 



O 




228° 



226° 224° 222° 

Galactic Longitude 



220° 



Fig. 7. A zoom of Fig. [s] showing the Canis Major dark 
complex. Toggle labels |. 



produced the final 2MASS/Nicer extinction map, shown in 
Fig. [3j To obtain these results, we smoothed the individual 
extinctions measured for each star, jA]^ |, using a moving 
weighted average 



(n) 



AK{e) 






(2) 



where Ak{0) is the extinction at the angular position 
and M^(") {0) is the weight for the n-th star for the pixel at 
the location 0. This weight, in the standard Nicer algo- 
rithm, is a combination of a smoothing, window function 
W{6 — 6^"-^^, i.e. a function of the angular distance between 
the star and the point where the extinction has to be in- 
terpolated, and the inverse of the inferred variance on the 
estimate of Ak from the star: 



M^(«)(6>) = 



T4^(6»-6»(")) 
Var(4")) 



(3) 



Note that the way Eq. (pi) is written, only relative values of 
W^^\ and thus of W[6 — 0^")) are important. Therefore, 
we can assume without loss of generality that the window 
function is normalized to unity according to the equation 



w{e)A^e = 1 



(4) 



For this paper, the smoothing window function W was 
taken to be a Gaussian with FWHM = 3 arcmin. Finally, 
the map described by Eq. ([2]) was sampled at 1.5 arcmin 
(corresponding to a Nyquist frequency for the chosen win- 
dow function). 

Similarly to Paper III, we also constructed a Nicest 
extinction map, obtained by using the modified estima- 
tor described in ,Lombardij (j2009_) . The Nicest map dif- 
fers significantly from the Nicer map only in the high 
column-density regions, where the substructures present in 
the molecular cloud produce a possibly significant bias on 
the standard estimate of the column density (see below). 
The largest extinction was measured close to LDN 1641 S, 
where Ak — 2.0 mag for Nicer and Ak — 5.0 mag for 
Nicest. 



Figures |4]-[7] show in greater detail the absorption maps 
we obtained for the Orion and Mon R2 star forming re- 
gions, the A Orionis bubble, the Rosette nebula, and the 
Canis Major complex. These maps allow us to better ap- 
preciate the details that we can obtain by applying the 
Nicer method to the 2MASS data. In the figures we also 
display the boundaries that we use throughout this paper 
and that we associate with the various clouds considered 
here. In particular, we defined 



Orion A: 
Orion B: 
Mon R2: 
A Orionis: 

Rosette: 
Canis Major: 



203° < I < 217° 
201° < I < 210° 
210° < I < 222° 
188° < I < 201° 
205° < I < 209° 
220° < I < 228° 



-21° <b< -17° , 
-17° < fo < -5° , 
-17° <fo< -7° , 
-18° <fo< -7° , 
-4° < fe < 0° , 
-4° < fe < 0° . (5) 



We note that the boundaries defined are somewhat arbi- 



trary, but as described below (see Sect. 3.2 and in particular 
Fig. [9]) there are strong indications that many (if not all) of 
the features within each of the defined regions are located at 
similar distances. We stress in any case that some of the re- 
sults presented (for example, the mass estimates discussed 
in Sect, m depend significantly on the chosen boundaries. 
The expected error on the measured extinction, a^ , 
(not shown here), was evaluated from a standard error 
propagation in Eq. (pi) [see below Eq. (pi)]. The main factor 
affecting the error is the local density of stars, and thus 
there is a significant change in the statistical error along 
galactic latitude. Other variations can be associated with 
bright stars (which are masked out in the 2MASS release), 
bright galaxies, and to the cloud itself. The median error 
per pixel for the field shown in Fig. [4] is 0.029 mag, while 
a significantly lower error of 0.019 mag is observed for the 
Rosette or Canis Major nebula. 

3. Statistical analysis 

3.1. Reddening law 

As shown in Paper 0, the use of multiple colors in the Nicer 
algorithm significantly improves the signal-to-noise ratio of 
the final extinction maps, but also provides a simple, direct 
way to verify the reddening law. A good way to perform this 
check is to divide all stars with reliable measurements in 
all bands into different bins corresponding to the individual 
A^^' measurements, and to evaluate the average NIR colors 
of the stars in the same bin. At first, this procedure might 
be regarded as a circular argument: in order to measure 
the extinction law one bins the colors of stars according 
to the estimated extinction! In reality, the method is well 
behaved and has interesting properties (Ascenso et al. 2011, 
in preparation): 

— The assumed reddening law is only used to bin the data, 
and is then iteratively replaced by the new reddening 
law as determined from the data. 

— This iterative process converges quickly and is essen- 
tially unbiased; the final reddening law essentially does 
not depend on the initial assumed one. 

— Simulations show that more standard techniques (such 
as binning in one simple color, e.g. J — K) suffer from 
significant biases mostly as a result of the heteroskedas- 
ticity of the data; additional complications are due to 
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Fig. 8. The reddening law as measured on various complexes. The plots show the color excess on J — i7 as a function of 
the color excess on H — K. Error bars are uncertainties evaluated from the photometric errors of the 2MASS catalog and 
taking into account also the intrinsic scatter in the i ndividual colors of the s tars, as derived from the control field. The 
solid line shows the normal infrared reddening law ( Indebetouw et al. 2005 1 . A significant discrepancy is observed for 
Mon R2 and possibly for A Orionis, while other clouds are in very good agreement with the normal infrared reddening 
law. 



the correlation of errors in the two colors J 
H- K. 



H and 



Figures [8] summarize the results obtained for a bin size 
of 0.02 mag in the various clouds. As shown by these plots, 
in al l clouds we have a good agreement between the stan- 
dard Indebetouw et al. ( 2005 ) infrared reddening law in the 
2MASS photometric system and the observed one. All plots 
show a systematic divergence below the reddening line at 
"negative" extincti ons, an effect due to the intr insic colors 
of dwarf stars (see Koornneef 1983 Bessell fc Brett 1988] ) . 
Specifically, the track of dwarf stars in the J — H vs. H — K 
color-color plot is slightly steeper than the reddening vec- 
tor (cf. Fig. [2]), and as a result stars in the bottom-left of 
the dwarf track, which in the simple Nicer scheme are as- 
sociated to negative extinction, exhibit an excess in their 
H ~ K color. This effect is not visible anymore as we go 
to positive extinction because the colors of these stars are 
then averaged together with the rest of the dwarf and giant 
sequence, and as imposed by the control field no systematic 
color excesses is present. 

Another common feature to all plots of Fig. |8]is an in- 
crease in the spread of points in the upper-right corner, i.e. 
for high column densities. As shown by the correspondingly 
larger error bars, this is expected and is a result of the low 
number of stars showing large extinction values (we recall 
that a constant bin size of 0.02 mag has been used every- 
where). Additionally, for some clouds, and in particular for 
Mon R2, we obse rve that the points fall sy stematically be- 
low the standard Indebetouw et al. (2005) reddening law. 
We believe that this difference is not due to intrinsic dif- 



ferences in the physical properties of the dust, but rather 
that it is caused by contamination from a population of 
relatively blue, young stars present in the field, both in 
the form of clusters and of dispersed OB associations (in 
particular, the area is known to host the Mon OBI associ- 
ation, see|Racine||1968[ [Herbst fc Racine||1976[ [Carpenter] 
fc Hodapp||2008r 

Finally, we stress that the check performed with Fig. |8] 
is a relative one: we can only verify that slope of the red- 
dening vector agrees wit h expectations from the standard 
Indebetouw et al. (2005) reddening law, but can not con- 
strain the length of the reddening vector. 

3.2. Foreground star contamination and distances 

Foreground stars are an annoyance for color extinction 
studies, since they dilute the signal coming from back- 
ground stars and add a source of noise to the maps. For 
nearby clouds, foreground stars are usually a small fraction 
of the total number of stars in all regions except in the 
very dense cores (where the density of background stars 
decreases significantly because of the dust extinction). In 
these conditions, one can safely ignore the bias introduced 
by foreground stars, and if necessary use Nicest measure- 
ments to alleviate the foreground bias in the high column- 
density regions. 

Foreground stars are easily selected in the dense regions 
of a molecular cloud as objects showing no or very little sign 
of extinction. This allows us to estimate with good accu- 
racy the local density of foreground stars in all high-column 
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Fig. 9. The local density of foreground stars, averaged on connected regions with extinction Af^^ > 0.6 mag. The three 
regions defined in Fig. |4] are shown again here. 
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Fig. 10. The distances of the clouds deduced from the density of foregrou nd stars. Solid lines show the density of 
foreground of stars as a function of the cloud distance, as predicted from the Robin et al. (2003) Galactic model. The 
grey areas correspond to the 95% (two-sigma) confidence regions for the density of foreground stars and the deduced 
confidence regions for the distances. 
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density clouds in our field. For this purpose, we selected 
connected regions with Ak{0) > 0.6 mag in our extinction 
map, and we computed for each of these the local den- 
sity of foreground stars (foreground objects were selected 
as stars with K-band extinction smaller than 0.3 mag, i.e. 
stars compatible with no or negligible extinction). Figure^ 
shows the results obtained in Orion and Mon R2; the av- 
erage numerical values obtained in the various regions con- 
sidered in this paper are also reported in Table [T] Two sim- 
ple results are immediately visible from Fig. [9 First, the 
Orion cloud shows significantly fewer foreground stars than 
the Mon R2 complex, as expected from the much larger dis- 
tance assigned in the literature to Mon R2. Second, all com- 
plexes have relatively uniform densities of foreground stars 
in their regions except the area around the Orion Nebula 
Cluster (ONC), marked with a circle in Fig. [9J This re- 
sult is also expected, since it is well known that the ONC 
region contains a high surface density of young stellar ob- 
jects (YSOs) embedded in the cloud that will contaminate 
our extinction measurements and a part of which will be 
considered as foreground stars by our selection criteria. In 
order to obtain unbiased results, we excluded this region in 
the analysis of the foreground star density. 

We note that a comparison of the densities reported 
in Table [l] with the average density of stars in the field, 
E ~ 9 000 deg^^, shows that the expected fraction / of 
foreground stars in the outskirts of the Orion complex is 
expected to be as low as / ~ 0.04, while a significantly 
larger value, / ~ 0.16, is expected for Mon R2. 

Interestingly, as shown in Paper III, we can also use 
foreground stars to estimate the distance of dark clouds. 
The technique relies on a comparison between the estimated 
density of foreground stars and the predictions of galactic 
models for the photometric dept h of the 2MASS ca talog. 
We used the Galactic model by Robin et al. (2003), and 
computed at the location of each cloud the expected num- 
ber of stars within the 2MASS pho tometric limits observed 
at various distances (see Fig. ]To|. The results obtained, 
shown in the last column of Table 111 presents several in- 
teresting aspects. First, the estimates for the Orion A dis- 
tance, when excluding the highly contaminated area of the 
ONC, is (371 ± 10) pc. This value is considerably less than 
the "standard" Orio n distance, 450 pc, but as discussed by 
(|2008|), this often quoted distance is actually 



Muench et al. 



reported by Genzel & Stutzki ( 1989 1 but it is not itself the 



result of any direct measurement, but more likely an aver- 
age of two different measurements. In general, in spite of 
the efforts over several decades to measure the distance of 



Complex 


iVf, 


/ 


Area 


Sfg 


Distance 








deg^ 


deg~2 


pc 


Orion A 


420 


0.047 


1.504 


279 ± 13 


371 ± 10 


Orion B 


298 


0.030 


0.936 


318 ±18 


398 ± 12 


A Orionis 


17 


0.018 


0.078 


219 ±53 


445 ± 50 


Mon R2 


207 


0.157 


0.144 


1430 ± 99 


905 ± 37 


Rosette 


275 


0.190 


0.083 


3330 ± 200 


1330 ± 48 


Canis Major 


99 


0.228 


0.036 


2730 ± 270 


1150 ±64 



Table 1. The average value of foreground stars A^fg and 
their fraction / found in the various complexes. Also re- 
ported in the last co lumn the estimated distance from a 
comparison with the Robin et al. (20031 Galactic model 
(see Fig.fTol). 



this cloud comple x, there is still a spre ad in recent estimates 
at the 10% level ( jMuench et al.||2008[ ). We stress, however, 
that our measurement is in good agreemen t with a recen t 



2007 



VLBI determination, (389 ±23) pc (Sandstr om et al. 
but lower than a measurement by |Menten et al.| 
(414 ±7) pc. 

Our Orion A distance is also very close to the Orion B 
distance, a result that is not unexpected given the physical 
relationship between the two clouds. 

We did attempt a measurement of the A Orionis dis- 
tance, in spite of the relative lack of dense material in the 
area. This paucity of high extinction material prevents us 
from estimating the density of foreground stars with the 
necessary confidence. Nonetheless, the result obtained is 
virtually identical to th e standard distance estimated by 
|Dolan & Mathieu (2001) using Stroomgren photometry of 
the OB stars in a main-sequence fitting in a theoretical H-R 
diagram, 450±50 pc. However, we are unfortunately unable 
to improve this result in terms of accuracy. 

The Mon R2 cloud is found to be at a significantly larger 
distance, (905 ± 37) pc, a result that compares well with 
the generally accepted dis tance of t his cloud, (830 ± 50) pc 



(Herbst & Racine 1976 



see also [Carpenter 
evident 



& Hodapp 



2008p . We stress that, as evident from Fig. M different 
clouds present in the MonR2 region seem to nave com- 
patible surface densities of foreground stars and there- 
fore distances, with the possible exception of the filaments 
close to NGC 2149. In order to test quantitatively this 
statement, we also considered subregions in the MonR2 
cloud. In particular, a measurement of the distance of the 
LBN 1015, LDN 1652, and LBN 1017 clouds (the "cross- 
bones") provides (750 ± 100) pc, while a measurement of 
the MonR2 core alone (i.e. the area around NGC 2170) 
gives (925 ± 150) pc. These two regions therefore seem to 
be at comparable d istances, a result tha t is in contradiction 
with what found by Wilson et al. ( 2005 1 . Unfortunately, the 



data we have are not sufficient to constrain the distance of 
the filaments around NGC 2149, and therefore we cannot 
securely associate them with MonR2 (nor with Orion A). 

The distance of the Rosette complex is generally 
obtained indirectly by measuring the distance of the 
NGC 2244 cluster. Typical values obtained from photomet- 



ric studies are around 16 50 pc (Johnson 1962 Perez et al. 
1 19871 [Park fc Sung||2002| , with the except ion of lOgura fc 
CafudaT'l^STl) that reports 1420 pc (see also Roman- Ziiiiiga 



& Lada 20 08| fo r a di scussion of these measurementsJT^ 
Hensbergc et ~al.| ( 2000 ) performed a spectroscopic analy- 
sis of the binary member V578 Mon, obtaining 1390 pc, a 
value that compares very well with our determination. 

Recent distance estimates of the Canis Major region 
(and in par ticular of its O B association) are close to 
1000 pc (see Gregorio-Hetem|[2b08 and references therein). 
Although our measurement appears to be slightly higher, 
the uncertainties of the data present in the literature and 
of our own determination are relatively large. 

In summary, our analysis shows that for clouds out- 
side the reach of Hipparcos parallaxes, a model-dependent 
distance obtained through number counts of foreground 
stars is a reasonable alternative of distance determination. 
Indeed, for all clouds studied here we obtained results that 
are in very good agreement with the data present in the 
literature (which often are obtained using dedicated obser- 
vations) . Additionally, the uncertainty we have is generally 
comparable to the one that can be obtained with compel- 
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itive techniques (with the notable exception of VLBI par- 
allaxes, which however are very expensive in terms of tele- 
scope time and need to rely on a secure physical link be- 
tween the source measured and the molecular cloud). The 
statistical error associated with this technique is linked to 
iVfg, the number of foreground stars observed in projection 
to the cloud. This quantity, for nearby clouds, increases lin- 
early with the cloud distance d, and therefore the relative 
statistical error on the estimated distance goes as l/\/d. 

3.3. Column density probability distribution 

Many theoretical studies have suggested that the turbulent 
supersonic motions that are believed to characterize the 
molecular clouds on large scales induce a log-nor mal prob 
ability distribution for the volume density (e.g. Vazquez 
T9941 [Padoan et al||1997b| [Passot fc Vazquez- 



Semadeni 



Semadeni 



T998| Scalo et al. 1998). Although this result 



strictly applies for the volume density, under certain as 
sumptions, verified in relatively "thin" molecular clouds, 
the probability distribution for the column density, i.e. 
the volume density integrated along the line o f sight, is 
also expected to follow a log- normal dist ribution ( Vazquez- 
Semadeni fc GarciapOOTI . Additionally, |Tassis et al.| ( |2010p 
recently have shown that log-normal distributions are also 
expected under completely different physical conditions 
(also plausible for molecular clouds) , such as radially strati- 
fied density distributions dominated by gravity and thermal 
pressure, or by a gravitationally-driven ambipolar diffusion. 
The log-normality of the column densities is verified 
with good approximation at low-column densities in many 
of the clouds that we have studied in the past and by simi 
lar res ults obtain ed by " 



Kainulainen et al.'(2009 



009]), 



Goodman 
et al.| ( |2009| ) , and lFroeEnch fc Rowlcs_ (2010| ^s shown by 
the plots in Fig. |11[ the region studied here confirms this 
general trend. The probability distributions of column den- 
sities for the various clouds were fitted with a log-normal 
distributions of the forrro 



hiA 



K) 



A 



K 



Ao 



exp 



(in(AK-Ao)-lnAi)' 
2a2 



(6) 



For some of the clouds, such as Orion B, A Orionis, and 
Mon R2, the fits appear to be better than for other ones, 
such as Orion A or Rosette. However, in all cases residu- 
als are well above the expected level^ and show system- 
atic and structured deviations even at low column densi- 



Note that the functional form used here differs, in the defini- 
tion of a, with respect to the form used in the previous papers. 
^ The theoretical error follows a Poisson distribution, and is 
therefore different for each cloud and each bin. In the range 



Table 2. The best-fit parameters of the four Gaussian func- 
tions used to fit th e column density probability distribution 
shown in Figs. 11 (see Eq. (rol) for the meaning of the various 
quantities) . 



Cloud 



Offset j4o Scale A\ Dispersion a 



Orion A 


-0.059 


0.193 


0.491 


Orion B 


-0.060 


0.145 


0.482 


Mon R2 


-0.163 


0.238 


0.214 


A Orionis 


-0.053 


0.121 


0.378 


Rosette 


-1-0.027 


0.079 


1.167 


Canis Major 


-1-0.013 


0.053 


0.978 
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Fig. 12. The A^ map of Eq. (l7| on the same region shown 
in Fig. [4] Note the significant increase observed in A^ close 
to the central parts of this cloud. 



ties. Additionally, all clouds show a positive residual at the 
higher column densities, approximately for A^ > 0.2 mag. 
The significance of these results and the goodness of the 
fits need to be further investigated. 

One perhaps surprising feature of Fig. [TT] is the pres- 
ence of a significant number of column density estimates 
with negative values. This could be either due to a zero- 
point offset in the control field or to uncertainties in the 
column density measurements, which naturally broadens 
the intrinsic distribution and possibly adds a fraction of 
negative measurements. Note also that the amount of neg- 
ative pixels observed is compatible with the typical error 
on our extinction maps, which is of the order of 0.03 mag. 

3.4. Small-scale inhomogeneities 



Lada et al. (1994) first recognized that the local dispersion 
of extinction measurements increases with the column den- 
sity. In other words, within a single "pixel element," the 
scatter of the individual stellar column density estimates 
is proportional to the average local column density esti- 
mate. This results implies the presence of substructures on 
scales smaller than the resolution of the extinction maps, 
and shows that theses substructures are more evident in 
regions with high column density. Substructures could be 
due either to unresolved gradien ts or to random fluctua- 
tions induced by turbulence (see Lada et al. 1999 ) 

The presence of undetected inhomogeneities is impor- 
tant for two reasons: (i) th ey might contain signatures of 
turbulent m otions (see, e.g. Miesch & Bally 1994 Padoan 
et al.|199"7a ), and (ii) they are bound to bias the extmction 
measurements towards lower extinctions in high-column 
de nsity regions (an d, especially, in the very dense cores; 



see 



Lombardi|[2009l ). 



In the previous papers of this series we have considered 
a quantity that traces well the inhomogeneities: 



A2(6>) 






(0) 



,(0)-(Var«O)W 



(7) 



displayed in Fig. |ll| the median error is approximately 0.1 mag, 
but since different bins are expected to be uncorrelated, the 
systematic offsets shown by the various clouds for Ak > 0.2 mag 
are highly significant. 
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Fig. 11. The probability distributions of pixel extinctions for the six cloud complexes. In each plot, the red, solid curve 
represents the best-fit with a log- normal distribution. Lower panels show the residuals with respect to the best-fit. 
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Fig. 13. The distribution of the A^ map as a function of 
the local extinction Ak for the map in Fig. [3j in logarith- 
mic grey scale. The dashed line shows the average values of 
A^ in bins of 0.05 mag in Ak- Note the rapid increase of 
A^ for Ak > 0.7 mag. As a comparison, the average vari- 
ance Var(y4.]^ ) on the estimate of Ak from a single star is 
approximately 0.033 mag^. 



The A^ map is defined in terms of the observed variance of 
column density estimates. 



.w 



EtiWi^ 



'[a'k' 



Me)]' 



e;=iW^(") 



(8) 



the average expected scatter due to the photometric errors 
and the intrinsic dispersion in the colors of the stars 



-1 



(0)^ 



_E„=i[W^*"He)] Var(i 



K ) 



E:=iW(-He) 



(9) 



and of the weighted average expected variance for the col- 
umn density measurements around 9 



(Var(4")))(0) 



E„M/(")Var(A^"0 



(10) 



As shown in Paper II, the combination of observables that 
enters the definition (l7| ensures that the expected value for 
A^ (9) is a weighted average of the square of local inhomo- 
geneities: 



(A^W) 






(11) 



The term inside brackets in the numerator of this equation 
represents the local scatter of the column density at 0^"' 
with respect to the weighted average column density in the 
patch of the sky considered Ak = E„ W^^"''^i^(^^"') / 
Sn W^''"''- Note also that this definition for Ak implies that 
the aver age of local inhomogeneities vanishes, as expected 
[cf. Eq. ml]: 










(12) 



Similar to the other papers of this series, we evaluated 
the A^ map for the whole field and identified regions with 
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large small-scale inhomogeneities. The A^ map for the re- 
gion considered in Fig. |4] is shown in Fig. [12] As usual, 
inhomogeneities are mostly present in high column den- 
sity regions, while at low extinctions (approximately be- 
low Ak < 0.4 mag) substructures are either on scales large 
enough to be detected at our resolution (2.5 arcmin), or are 
negligible. 

The average A^ as a function of the local extinction, 
Ak: for the map in Fig. l6]is presented in Fig. 13 A com- 
parison of the dashed line, representing the average value 
of A^ in bins of 0.05 mag in Ak, with the average variance 

Var(v4]^ ) on the estimate of Ak from a single star, which 
is approximately 0.033 mag^, shows that unresolved sub- 
structures start to be the prevalent source of errors in ex- 
tinction maps for Ak > 0.5 mag. As discussed in |Lombardi| 
( 2009 1, this column density value gives also an approximate 



upper limit for the robustness of Nice and Nicer extinc- 
tion studies, since both these methods are based on the 
implicit assumption that the extinction is uniform within a 
resolution clement (in our case, within the region where the 
window function W{0) is significantly different from zero). 
In contrast, the Nicest algorithm is explicitly designed to 
cope with unresolved substructures, and therefore it should 
be used in the region considered when Ak > 0.5 mag. We 
stress that the use of Nicest does not remove the unre- 
solved inhomogeneities (and thus does not make the A^ 
map flat), but rather it makes sure that the estimate of ^i^- 
is unbiased even if A^ is non- vanishing. 

4. Mass estimate 

Masses of the clouds were evaluated using the standard 
relation 



M^d^HpK / AK{e)d^e 



(13) 



where d is the cloud distance, fi is the mean molecular 
weight corrected for the helium abundance, and Pk — 
1.67 xlO^'^ cm-^ mag-i is the ra tio \N{m)+2N{B.2)\IAK 
JSavage fc Mathis|[T979| see also |Lilley|[l955| jBohliii et al. 
i978p . Assuming a standard cloud coinposition (63% hy- 
drogen, 36% helium, and 1% dust), we find /x = 1.37. The 
integral above was evaluated either over the whole area of 
each cloud, as defined in Eqs. ([5|, or inside contours above 
a given extinction threshold. This latter option allowed us 
to avoid using the "total" mass of a cloud, which is ill de- 
fined and strongly depends on the contours chosen, in favor 
of the mass of a cloud inside specified extinction contours. 



Cloud 



Distance Mass (Mq) 

Tbtal Ak > 0.1 Ak > 0.2 



Orion A 
Orion B 

Mon R2 
A Ori 
Rosette 
Canis Major 



371 pc 
398 pc 
905 pc 
445 pc 
1330 pc 
1150 pc 



75 700 
95100 
392 000 
102 000 
233 000 
205 000 



66 400 

68 300 

222 000 

47400 

193 000 

137000 



45100 
36100 
73 300 
11500 
137000 
76 900 



Table 3. The masses of all complexes studied in this pa- 
per. The table reports the assumed distance (cf. Table I]), 
followed by the estimated masses over the whole area, and 
over areas within the Ak > 0.1 mag and the Ak > 0.2 mag 
contours. 




Fig. 14. The cumulative mass enclosed in iso-extinction 
contours for the various clouds. All plots have been ob- 
tained from the extinction map shown in Fig. l6J and 
have thus the same resolution angular limit {FWHM = 
3 arcmin). 
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Fig. 15. Same as Fig. 14 but the extinction maps for the 
various clouds were degraded to the same physical resolu- 
tion FWHM -1.13 pc. 



We stress that the masses indicated for the various 
clouds refers to the boundaries taken in this paper. In par- 
ticular, if we increase the southern boundary of MonR2 
from b > —17° to b > —15°, thus avoiding the filament 
south of NGC 2149, then the mass estimates for MonR2 de- 
crease to 311 000 M0 (total), 166 000 M© {Ak > 0.1 mag), 
and 58 600 Mq {Ak > 0.2 mag). Correspondingly, if this 
area encloses Orion A, then the masses of this complex in- 
crease to 86 000 Mq (total), 75 800 Mq {Ak > 0.1 mag), 
and 47 600 Mq {Ak > 0.2 mag). 

Uncertainties on the cloud masses are difficult to evalu- 
ate, because they are essentially dominated by systematic 
errors. Note that the errors on the mass due to statisti- 
cal noise of the extinction maps are negligible. Instead, the 
major sources of errors are the uncertainty on the cloud 
distances (which, because of the dependence of the mass 
on the square of the distance, can be relevant), and pos- 
sible systematic errors on the zero-point of the extinction 
maps (typically due to a residual extinction present in the 
control field). 
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As in the other papers in this series, it is interesting 
to investigate the cumulative mass as a function of the ex- 
tinction threshold. This plot, shown in Fig. [14J provides a 
simple measure of the structure of the molecular clouds, and 
of the relative importance of low- and high-density regions 
within each cloud. Additionally, as shown by |Lada et "ah 
(2010), the cloud mass at high densities (in particular, for 
Ak > 0.8 mag) appears to correlate with the overall star 
formation rate of the cloud. 

Figure [15] shows the same plot of Fig. [Mj but using 
this time the same physical resolution for all clouds. That 
was accomplished by degrading the extinction map of the 
nearby clouds to match the physical resolution of the most 
distant cloud, the Rosette. This procedure makes even more 
evident the difference between the various complexes; inter- 
estingly, the smoothing applied also strengthens differences 
between clouds at the same distance, such as Orion A and A 
Orionis. This last point highlights differences in the density 
structure present in the clouds studied here. In particular, 
clouds that present a plateau of relatively high values of 
extinction are not affected too much by the smoothing ap- 
plied in Fig. [15] on the contrary, clouds such as A Orionis 
that have dispersed material will have their mass dispersed 
over an even larger area, and will therefore present sharp 
decreases in their cumulative mass function for relatively 
small values of Ak ■ 



5. Conclusions 

The following items summarize the main results presented 
in this paper: 

— We measured the extinction over an area of ^ 

2 200 square degrees that encompasses the Orion, the 
Monoceros R2, the Rosette, and the Canis Major molec- 
ular clouds. The extinction map, obtained with the 
Nicer and Nicest algorithms, has a resolution of 

3 arcmin and an average l-cr detection level of 0.3 visual 
magnitudes. 

— We measured the reddening laws for the various 
clouds and showed that they agree with the standard 



Indebetouw et al. (2005) reddening law, with the ex- 
ception of the Mon R2 region. We interpret the discrep- 
ancies observed as a result of contamination from a pop- 
ulation of young blue stars (likely including members of 
the Orion OBI association). 

We estimated the distances of the various clouds by 
comparing the density of for egroun d stars with the pre- 
diction of the Robin et al. (|2003[) Galactic model. All 



values obtained were found to be in very good agree 
ment with independent measurements. 
We considered the column density probability distribu- 
tions for the clouds and obtained reasonable log-normal 
fits for all of them. 

We measured the masses of the clouds and their cumu- 
lative mass distributions and found differences in the 
internal structure among the clouds. 
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